On the QCD phase diagram at finite chemical potential 
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Abstract. We present results for the phase diagram of QCD with two massless quark flavours as obtained from a first-prin- 
ciples functional renormalisation group approach. In particular we compute order parameters for chiral symmetry breaking 
and quark confinement at vanishing and non-zero imaginary chemical potential. Our analytical and numerical results suggest 
a close relation between the chiral and the deconfinement phase transition. We discuss the properties of dual order parameters 
at imaginary and real chemical potential. 
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One of the unresolved problems in high energy 
physics is the structure of the QCD phase diagram. This 
concerns in particular the transition(s) between a decon- 
fined and chirally symmetric phase with microscopic de- 
grees of freedom, quarks and gluons, to a phase of col- 
our-neutral macroscopic bound states, hadrons, with bro- 
ken chiral symmetry. The deconfinement transition is re- 
lated to the breaking of the center symmetry of the gauge 
group and is driven by gluon dynamics. The chiral tran- 
sition is triggered by strong gluon-induced quark inter- 
actions. We report on results of the QCD phase diagram 
and the relation between the two phase transitions [1]. 
The question whether both transitions are related is sub- 
ject of an ongoing debate, see e.g. Refs. [2, 3] for lattice 
and [4, 5, 6, 7, 8, 9] for model studies. 

Approach. We compute the QCD effective action 
with fRG techniques, which allow us to include all quan- 
tum fluctuations step- wise at each momentum scale, for a 
recent overview see [10]. This is achieved by integrations 
over small momentum shells, generating a flow from 
the microscopic action in the UV towards the macro- 
scopic action in the IR. In our two-flavour calculation in 
the chiral limit [1] we include the Yang-Mills sector of 
QCD [11, 12] and the matter sector [13, 14, 15, 16] and 
couple them via dynamic quark-gluon interactions. This 
approach has already been applied to the chiral phase 
boundary in one-flavour QCD at finite chemical poten- 
tial [13]. The confining properties are included via the 
full momentum dependence of the ghost and gluon prop- 
agators [11, 12, 17] and the matter sector incorporates 
dynamical mesonic degrees of freedom. 

Order parameters. The order parameter for the de- 
confinement phase transition is the Polyakov loop; it is 
proportional to the energy needed to put a quark into the 



theory. In our study we implement the Polyakov loop as 
defined in Ref. [11, 18]. The order parameter for the chi- 
ral transition is related to the quark condensate. 

Recently so-called dual order parameters for the de- 
confinement phase transition have been defined [19]. 
This has been extended in [20] to any observable that 
transforms non-trivially under center transformations 
and has been applied in Refs. [21, 22, 23, 24, 25]. An 
element of the center Z(N C ) of the gauge group is given 
by z = \e 2%iQ \ where Z = 0, 1/3,2/3 for SU(3). It fol- 
lows immediately that the sum over all center elements is 
zero in the symmetric phase and non-zero in the broken 
phase. This means that any observable that transforms 
non-trivially under center transformations is an order pa- 
rameter for confinement. 

We have extended the above setting to imaginary 
chemical potential. It can be incorporated in generalised 
boundary conditions of the quarks and rewritten in terms 
of physical quarks with anti-periodic boundary condi- 
tions 

yr e (x) = with yr(x) = yr e= o(x) (1) 

and j3 = l/T. Due to the periodicity in the angle 0, 
we can Fourier decompose general observables G$ = 
(@[e 2Kiet IP if]) which depend on the quark fields 

0e= JV^O/. (2) 

zez 

Under center transformations the 0/ are multiplied with 
a center element, 0\ — >> z l Oi and hence are order param- 
eters for confinement as they are proportional to the sum 
over all center elements and thus vanish in the symmet- 
ric phase. One example is = e 2m0 (L) , where L is the 
Polyakov loop variable. In particular for / = 1 we find 

&= ( l dQe- 2nie @ e = C d0O x =O x . (3) 
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FIGURE 1. The pion decay constant as a function of imagi- 
nary chemical potential and temperature. 



Thus the dual Polyakov loop in QCD is (L) . 

In general, observables @q can either be evaluated in 
QCD with anti-periodic quarks, see e.g. Refs. [19, 20, 21, 
22, 23, 24, 25], or in QCD at imaginary chemical poten- 
tial, QCDe, with 6— dependent boundary conditions [1]. 

The Dirac action with quark fields defined in (1) reads 

J \jf e (0 + im)\ife, (4) 

where ij) — $ — igji. This can be rewritten such that we 
obtain an additional term which has the same form as an 
imaginary chemical potential 6 : 

J Y^ilfi + im-ln^ybO^ W- (5) 

Imaginary and real quark chemical potential are re- 
lated via 6 = —ijifi/lTZ. The effective action of QCD 
is then periodic under a transformation of 6 — ^ 6 + 6 Z \ 
QCDfl = QCD 0+ , as the transformation of 6 is can- 
celed by the center transformations of the fields. Observ- 
ables G$ related to the effective action show the same 
periodicity, namely the Roberge- Weiss (RW) periodic- 
ity. However, this also means that if RW periodicity is 
not broken explicitly, all 0\ vanish. The QCD phase di- 
agram at imaginary chemical potential shows a smooth 
transition until = 1/6 but then it displays a discontinu- 
ity: the Polyakov loop RW phase transition at Trw [26]. 

Presence of a fixed background field. The RW pe- 
riodicity of the generating functional is broken in the 
presence of a current /. The dual observables 0\ then 
no longer vanish and as derived above, serve as order pa- 
rameters for confinement. This is implemented by a 6- 
independent gauge-field background (p. One example of 




FIGURE 2. AP(r, 6) as a function of temperature and imag- 
inary chemical potential. 



o.o 




FIGURE 3. The dual quark mass parameter as a function of 
temperature and imaginary chemical potential. 



the resulting order parameters is the dual density, which 
is proportional to the logarithm of the generating func- 
tional. Therefore it grows like T 3 at high temperatures as 
it is proportional to the first moment of the grand canon- 
ical potential. Integration by parts yields the fermionic 
pressure difference AP(T, 0) = P(T, 0)-P(T,0). 

As the &q are observables in different theories, distin- 
guished by the boundary condition, the & vanish only if 
QCDfl is in the center symmetric phase for all boundary 
conditions. 

Fig. 2 displays the fermionic pressure difference in the 
presence of a 6 -independent gauge-field background (p. 
The RW symmetry is broken, however 0—^0 + 1 still 
holds. Instead of imaginary chemical potential, one can 
think of this as an order parameter for different theories 
which are distinguished by the value of the angle 0. 

Another dual observable is the dual quark mass param- 
eter M[(j>j], see Eq. (3) and Me[<j>j] in Fig. 3. Due to the 
presence of a fixed background field it does not vanish 
in the broken and in the symmetric phase. At = 1/2 
closely above the phase transition it increases with y/f 
and then linearly as the quarks effectively have bosonic 
Matsubara frequencies at 6 = 1/2. At = the dual 
quark mass is zero above the transition. For vanishing 
current it is related to the pion decay constant in QCD: 
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FIGURE 4. The pion decay constant, the dual density and 
the Polyakov loop as functions of temperature, %l = djL, 
Zdual = d T n. 
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FIGURE 5. The QCD phase diagram at imaginary chemical 
potential. The grey band represents the width of %L- Black dots 
indicate the endpoints of the Polyakov loop RW transitions. 



the slice of the 3D plot at 6 = is the (normalised) line 
shown in Fig. 4. 

Results. In our calculation we include the back-reac- 
tion of the matter sector on the gauge sector. Moreover, 
we do not use input from, e. g., lattice calculations to 
model the gauge dynamics. In other words, the gauge 
and the matter sector as well as their interplay are treated 
self-consistently within our approach. 

At vanishing chemical potential we consider the order 
parameters of the chiral and the deconfinement phase 
transition, the Polyakov loop, the dual density and the 
pion decay constant, see Fig. 4. 

Above T c — 180 MeV the pion decay constant van- 
ishes and chiral symmetry is restored. The dual den- 
sity and the Polyakov loop both show a peak in their 
temperature derivative at w 178 MeV. This provides a 
non-trivial consistency check of our approximation as 
the Polyakov loop is computed from gluonic correlation 
functions, whereas the dual density is computed from 
matter correlation functions. We find that the chiral and 
the deconfinement transition agree within a few MeV. 

Fig. 1 displays the pion decay constant as a function 
of imaginary chemical potential and temperature. For 
T > T c % it vanishes and it is non-zero below T c %. The 
RW symmetry is found. Moreover we find a second order 
phase transition as expected in the chiral limit. 

Fig. 5 is a plot of the QCD phase diagram at imag- 
inary chemical potential. The chiral and the deconfine- 
ment transition agree within the width of the temperature 
derivative of the Polyakov loop throughout the phase dia- 
gram. The deconfinement transition occurs at lower crit- 
ical temperatures than the chiral transition. This was also 
found by lattice computations [27, 28]. First results indi- 
cate that this persists at real chemical potential [29] . 
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